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ABSTRACT 

We present a catalog of 17 filamentary X-ray features located within a 68'x34' 
view centered on the Galactic Center region from images taken by Chandra. These 
features are described by their morphological and spectral properties. Many of the X- 
ray features have non-thermal spectra that are well fit by an absorbed power-law. Of 
the 17 features, we find 6 that have not been previously detected, 4 of which are outside 
the immediate 20'x20' area centered on the GC. 7 of the 17 identified filaments have 
morphological and spectral properties expected for pulsar wind nebulae with X-ray 
luminosities of 5x10'^^ to 10'^'* ergs s~^ in the 2.0-10.0 keV band and photon indexes 
in the range of r=l.l to 1.9. In one feature, we suggest the strong neutral Fe Ka 
emission line to be a possible indicator for past activity of Sgr A*. For G359. 942-0. 03, a 
particular filament of interest, we propose the model of a ram-pressure confined stellar 
wind bubble from a massive star to account for the morphology, spectral shape and 6.7 
keV He-like Fe emission detected. We also present a piecewise spectral analysis on two 
features of interest, GO. 13-0. 11 and G359.89-0.08, to further examine their physical 
interpretations. This analysis favors the PWN scenario for these features. 
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1 INTRODUCTION 

The Galactic Center (GC) region is a unique environment 
that is home to numerous energetic processes. Radio, in- 
frared and X-ray observations show structures that can only 
be distinguished in our own galaxy due to its proximity. The 
GC thus provides a laboratory in which to probe the Galac- 
tic nuclear environment and the interactions between star 
formation regions, the interstellar medium (ISM) and the 
supermassive black hole Sgr A*. These observations can aid 
in understanding the nuclear environments of nearby galax- 
ies. 

Radio observations have detected prominent non- 
thermal radio filamen ts that ha ve been intensively stud- 
ied (e.g.[Yusof-Zadch, Morris fc C hance 1984). Many of the 
non-thermal radio filaments may come from milliCauss mag- 
netic fields illuminated by energetic particles. These fields 
are seen to be primarily perpendicular to the galactic plane 
and indicate the general structure of the magnetic field 
around the GC. The exact origins and implications of the 
radio non-thermal filaments in relation to the GC are still 
in debate. 

Chandra observations have detected numerous large 
scale X-ray structures from thermal to non-thermal fil- 



aments (e.g iBamba et al. I l2002l : ILu. Yuan, fc Lou I l2008l : 
iMuno et al. Il2007h primarily in the local 20'x20' area cen- 
tered on the GC. The commonly accepted scenarios for non- 
thermal thread-like X-ray fea tures in the GC include pulsar 
wind nebulae (PWNe) (e.g. IWang. Li. fc Begelman I Il993l : 
IWang. Lu fc Gotthelf "1120061 '). with the filament-Uke struc- 
ture produced through ram-pressure confinement or strong 
magnetic fields, and magnetohydrodyn amical shock fronts 
from supernova remnants (SNRs) (e.g lYusef-Zadeh et al. I 
l2005h . For the SNR case, the elongation of features rep- 
resent shock fronts in the ISM. The average non-thermal 
spectrum of a typical PWN is generally well fit by a power 
law model with a photon i ndex in the range of r=l.l- 
2.4 (jCotthelf fc Olbert |[20o3 ') and X-ray luminosities in the 
range of 10^'^ to 10^^ ergs s~^ fo r the 0.2-10 keV energy band 
(e.g. ICaensler fc Slane I l2006l : H aspi. Roberts fc Harding I 
llOOtf). The PWN scenario is motivated primarily by the 
similarities between the features seen in the GC and 
X-ray features associated with known pulsars (see also 
iKargaltsev fc Pavlov 1 12008| ) . Few cases of Ka transitions 
from neutral iron (Fe) are proposed to stem from the re- 
flection of hard X-rays from extern al sources, such as th e 
supermassive black hole Sgr A* fe.g. iKovama et al. |[l989l ). 
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By analyzing these X-ray features, one can in principle trace 
gas dynamics and magne tic fields in the GC region (e.g. 
IWang. Lu . fc Lang ' '2002*1 and examine the history of GC 
activity (iKova ma et al. .,1989 ). 

Previous broad scale studies of X-ray features fo- 
cus on the immediate 20'x20' area surroundi ng Sgr A* 
(e.g. iMuno et al.1l2007l : ILu. Yuan, fc Loull2008h . Here, we 
present a wide scale study on 17 X-ray thread-like features. 
The features are detected over a roughly 68'x34' area of the 
sky, relating to approximately 158 pc by 79 pc, around Sgr 
A* based on all available observations. While some of these 
features have been examined in previous literature, the im- 
proved counting statistics of this study allows for piecewise 
analysis that will be able to detect trends in the spectra 
and therefore better co nstrain the physical model (see also 
ILu. Yuan, fc Loull2008l ). In addition to the piecewise anal- 
ysis of some features, we also present an additional physical 
model in order to explain a feature with a Helium (He)-like 
Fe Ka emission line, comet-like morphology and steep spec- 
trum. For consistency within this study, the spectral and 
morphological characteristics were performed for all features 
listed. All physical distance measurements assume an 8 kpc 
distance to the GC and all errors, unless otherwise stated, 
are given to the 90% confidence level. 



2 OBSERVATIONS AND DETECTION 

The Chandra X-ray Observatory observed the GC with 
the Advanced CCD Imaging Spectrometer imaging array 
(ACIS-I) for a total integrated exposure time of ^2 Msec 
betw een 1999 and 2007 over 81 observations (|Muno et al. I 
I2OO8I ). The array is composed of four individual CCDs 
that operate together to provide a 17'xl7' view with sub- 
arcsecond resolution at the center of the array. The images 
were first reprocessed using the standard Chandra Interac- 
tive Analysis of Observations CIAO routines (version 3.4.0) 
and then combined by celestial coordinates to produce a 
merged events file. 

Initial identification of X-ray filaments comes from a 
careful visual examination of the merged event image in the 
2.0-8.0 keV band. We examine the merged events file under a 
logarithmic scale with different upper/lower limits and look 
for extended, linear structures. From this initial identifica- 
tion, we removed features whose morphologies and apparent 
extended structures are not well confined with respect to the 
feature in question, leaving only those with tightly correlated 
extended X-ray emission. Within regions of strong diffuse 
X-ray emission, e.g. the Sgr A* complex, we limited our se- 
lection to well constrained extended X-ray emission regions 
with size scales larger than 3" x 1" . These morphologically 
selected features were culled once more after taking their 
individual spectra (see below for details on spectral extrac- 
tion) and rejecting them if their spectra appear consistent 
with their local background; e.g. in spectral shape. Based on 
our selection criteria, these features are thus X-ray bright, 
in respect to their local backgrounds, and have well defined 
linear morphologies; however, due to the complexity of X- 
ray emission around the GC and the uneven sensitivity of 
this survey, this sample is far from complete and represents 
only a fraction of the extended X-ray features indigenous to 
the GC. Fig. 1 shows the composite intensity map in the 2.0- 
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Figure 1. Full view of the 81 compiled observations taken by 
Chandra around the GC. The center square identifies the 20'x20' 
region of Fig. 2. Features outside this region are selected by their 
surrounding ellipses and labeled based on their galactic coordi- 
nates. The image is the exposure-corrected intensity map in the 
2.0-8.0 keV band produced by dividing the ACIS-I count map 
with the exposure map. This and subsequent images are given in 
logarithmic scaling with North defined as up and East as left. 



8.0 keV band with identified features labeled and indicated 
by their surrounding ellipses. Fig. 2 shows the inner 20'x20' 
around Sgr A* with features again labeled and constrained 
by their surrounding ellipses. In extracting their source and 
background spectra, point sources and other such possible 
contaminators were excluded from the count and exposure 
maps. We use local background subtraction for each fila- 
ment due to the widespread area and time integration over 
which the filaments are found. For those that are surrounded 
by strong diffuse emission, e.g. within the Sgr A* complex, 
backgrounds were selected to account for the immediate in- 
tensities. 

Typical spectral extraction methods pose a problem 
for the features that were covered by multiple observations. 
While some of the 81 observations have the same celestial 
coordinates in their pointing, they could have different roll 
angles. As such, a feature could fall into different positions 
on the ACIS-I detector over different observations including 
the gap between the CCDs that compose the array. Addi- 
tionally, a feature may be completely contained within one 
observation while partially covered in another. Therefore, 
the traditional spectral extraction could not be used due to 
time and spatial variability of the Auxiliary Response File 
(arf) and Redistribution Matrix File (rmf) of the instru- 
ment. Another way is to extract the spectra of the filaments 
from individual observations and use "ftools" to merge 
them into a final spectrum. However, this is just suitable for 
a single object and can become difficult when dealing with 
many filaments over many observations. Therefore, we de- 
veloped a new method to directly extract the spectrum from 
the merged events file. We first reproject all the events in 
the merged events file of a certain filament into the instru- 
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Figure 2. Magnified view of the intensity map from Fig.l on 
the inner 20'x20' area centered on Sgr A*. Identified features are 
again selected and labeled by their surrounding ellipses. 



ment coordinate, since the calculation of the ARF and rmf 
is based on events' positions on the CCD, not the celestial 
coordinate. The ARF is not only a function of the events' 
positions, but also depends on when the observation was 
taken. The quantum efficiency of the CCDs has been seen 
to degrade over time in part due to molecular contamination 
build up on the the CCD itself and/or the filter. However, 
we notice that above 2 keV, the ARF of ACIS-I decreases less 
than 5% during the 7 years observational period (absorption 
towards the GC is large for energies below 2 keV, making 
the 2 keV lower limit optimal in terms of spectral analysis) . 
Using this lower limit, the ARF can then be considered to be 
time-independent. The rmf file is very stable and is just a 
function of position, not time (we have checked the rmf file 
used by the CIAO routines for each observation one by one. 
All of them use the same calibration file). Therefore, a new 
response matrix file, rsp (the multiple of the ARF and rmf 
file), for each small grid on the CCD (the place of each grid 
in the CCD is predefined by the calibration file) with events 
was created from the calibration file and the total number 
of events falling into each grid was used as the weight for 
merging all the new RSP files in different grids into the final 
merged rsp. We have tested this method in the center 7' re- 
gion around SGR A* and found that the fitting result from 
this method is nearly identical to that from the same region 
in the observation with the longest exposure time, proving 
that our method is robust and reliable. 



3 ANALYSIS OF X-RAY FEATURES 

For each of the identified filaments, we classify its basic ap- 
pearance as well as spectral proprieties. We broadly classify 
the features into two morphological categories: 'filamentary' 
and 'cometary'. A feature is considered to be cometary if 
its morphology appears to taper off in one direction, sim- 
ilar to a comet's tail, or if it has apparent asymmetry in 
its brightness across its major axis. Features that do not 
meet these criteria are defined as filamentary, this includes 
features that appear to be brightest in the center. This mor- 
phological classification along with the sizes of the features 



are approximations based on the observed flux and appar- 
ent surface brightness relative to the local background. The 
name associated with each feature comes from the approxi- 
mate center or brightest point of the feature. 

Fig. 3 through 7 show the individual features with el- 
lipses indicating regions used in extracting the source and 
background spectra. Linear X-ray features can be seen in 
some images in addition to the features identified. Many of 
these were either not selected due to the possibility of be- 
ing small scale structures of a larger X-ray emitting region 
or rejected based on our selection criteria. The plot next 
to each Chandra count map gives the extracted spectrum 
for each feature. In each case, the spectrum is well fit us- 
ing a power law model with foreground absorption in xspec 
(pha(po) for short). A Gaussian model is included for those 
features which exhibit line emission of 6.4 keV neutral Fe or 
6.7 keV He-like Fe (pha(po+GA)). The spectral parameters 
for the pha(po) or pha(po-I-GA) single fit models as weU 
as the morphological properties for the features are listed in 
Tables 1 and 2. It is well known that extinction towards the 
GC can vary g reatly over scales on order of arc-seconds (e.g. 
IScoville et al. |[2003i l. As such, each feature is modeled with 
an independent value of column density (A'h) rather than 
a single fixed A^h for all features based on the assumed 8 
kpc distance to the GC. The following summary of features 
focuses on those with unique traits or observations. 



3.1 GO. 223-0. 012 

GO. 223-0. 012 (Fig. 3a) is located near the bright X-ray point 
source CXOGCS J174621.05-284343.2, removed from the 
image in Fig. 3a with its location marked, to the northeast of 
the feature. The feature's position as well as orientation re- 
jects the possibility that GO. 223-0. 012 could be the result of 
a CCD readout streak ft-om CXOGCS J174621.05-284343.2. 



3.2 GO.13-0.11 

GO.13-0.11 (Fig. 3b) was previously studied by 
IWang. Lu. fc Lang I (|2002l ') as a possible PWN. Unlike 
other PWN candidates in the GC region, G0.13-0.11's 
morphological proprieties are fairly unique in that it is 
curved analogous to that of a co-existin g non-thermal radio 
filament in the Radio Arc region (see IWang. Lu, fc Lang I 
|2002| . fig. 1). The morphology of GO.13-0.11 appears to 
include two tail-like segments stemming from the apparent 
point source giving the feature a wing-like appearance. 
We identified the feature as cometary due to th e more 
prominent eastern segment. IWang. Lu. fc Lang I (|2002l ) 
proposed that the unique morphology of this feature is 
likely due to interactions of high-energy particles from 
the pulsar in a strong magnetic field, traced by the radio 
polarization of the Radio Arc region. To determine if there 
is some spectral evolution across the feature, as expected 
for PWNe, we perform a piecewise analysis of the apparent 
point source and the 'wings' of GO.13-0.11 (see Fig. 8a). In 
separating the point source from the wing, we adopt the 
1.5" extraction radius used in lWang. Lu. fc Lang I (|2002l ). A 
joint fit of the point source and total wing spectra, keeping 
A^'h a common parameter between both features, produces 
a good fit (xVd-o.f=79.7/86) with 7Vh = 
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Figure 3. Left-hand panel: the X-ray count maps of individual features in the 2.0-8.0 keV band. The upper and lower limits for each 
count map are 12 and 0.5 counts/arcsec^ for GO. 223-0. 012 (a), 23 and 2.3 counts/arcsec^ for GO. 13-0. 11 (b), 69 and 6.9 counts/arcsec^ 
for GO. 03-0.06 (c), and 207 and 4.1 counts/arcsec^ for GO.017-0.044 (d), respectively. Source regions are identified by the solid ellipses 
with background regions enclosed by dashed ellipses. In 3a, the empty region results from the removal of events around the bright X-ray 
point source CXOGCS J174621. 05-28 4343. 2 whose position is given by the small circular region. Fig. 3c and 3d include the 'West Clump' 
region from Yuan, &: Lou ] 1I2OO8I) , labeled as 'WC'. Right-hand panel: the X-ray spectra for the matching features given in the left 
column in the 2.0-8.0 keV band. All spectra are fit with an absorbed power-law model (pha(po)). The model for GO.017-0.044 is fitted 
with a Gaussian centered at 6.4 keV in addition to the pha(po) 



cm~^ and F of l.Ol g for the point source and I.IIqI 
for the wing; these val u es rou ghly agree with the analysis 
of IWang, Lu. fc Larigl l|2002l ). Segmenting the wing and 
performing the joint fit again provides a joint fit column 
density of NH=5.2tls x cm"^ and F of O-dtgi, 

O.bto'g, and 1.2to'6 for the point source, wing segment 
containing the point source and wing segment away from 



the point source, respectively, with x^/d.o.f=90.7/89. These 
best fit spectral parameters determined from the joint fits 
are in agreement with those listed in Table 2 which were 
determined for the feature as a whole. While the photon 
index appears to have steepened away from the point 
source, the error bars for the photon indexes of the two 
wing segments have considerable overlap. As such, we can 
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Figure 4. Continued X-ray count maps and spectra similar to Fig. 3. Upper and lower limits for each count map are 103 and 4.1 
counts/arcsec^ for GO.007-0.014 (a), 410 and 4.1 counts/arcsec^ for G359. 97-0.038 (b), 207 and 4.1 counts/arcsec^ for G359.97-0.009 (c), 
and 410 and 12 counts/arcsec-^ for G359. 964-0. 052 (d), respectively. 



not fully constrain the spectral evolution with the available 
data. 



3.3 GO. 017-0. 044 

Spectral analysis of GO. 017-0. 044 (Fig. 3d) shows a strong 

6.4 keV emission line of neutral Fe Ka transition. Given 
its projected separation from Sgr A*, approximately 11 pc, 
and this strong neutral Fe emission, possible explanations for 
this feature include the reflection of radiati ve illumination 
from Sgr A* past activity as proposed by iKovama et al. I 



(| 19891 ) and collisions in molecular clouds between low en- 
ergy cosmic-ray electrons and ions (e.g. lValinia et al~ll2000l ). 
though it is difficult to say with certainty due to the 
low photon counts. The relative flatness of the spectrum 
indicates that the feature may result from a high tem- 
perature plasma. This, however, is inconsistent with the 
weakly ionized Fe responsible for the 6.4 keV line unless 
the plasma is in a highly non- c ollisio nal ionization equi- 
librium state. ILu. Yuan, fc Lou~l l|2008h studied a clump of 
diffuse X-ray emission located to the west of GO. 03-0. 06 
(Fig. 3c) which has similar line emission and equivalent 
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Figure 5. Continued X-ray count maps and spectra. Upper and lower limits for each count map are 207 and 8.3 counts/arcscc^ for 
G359.96-0.028 (a), 1200 and 21 counts/arcsec^ for G359.95-0.04 (b), 207 and 12 counts/arcsec^ for G359.942-0.03 (c), and 410 and 4.1 
count/arcsec^ for G359.94-0.05 (d), respectively. G359. 942-0. 03 is fitted with a 6.7 keV Gaussian in addition to the pha(po) model. 



width as GO.017-0.044. This 'West Clump' (as referred to by 
ILu. Yuan, fc Loull2008t ) is seen in the images for GO. 03-0. 06 
and GO.017-0.044 (Fig. 3c and 3d, respectively). In addition 
to the 6.4 keV line similarities, the photon index of GO.017- 
0.044 is also comparable to that of the 'West Clump' albeit 
with fewer photon counts and thus larger 90% confidence 
level error bars. The similarities between the spectra of this 
West Clump and GO.017-0.044 indicate that they may be 
produced through the same mechanisms. 



3.4 G359. 942-0. 03 

G359. 942-0. 03 (Fig. 5c) displays morphological properties 
similar to a ram-pressure confined PWN but with a no- 
ticeable 6.7 keV emission line from He-like Fe. The pres- 
ence of this line and the steepness of the spectrum indi- 
cates that the emission is likely thermal in nature, contrary 
to the typical non-thermal PWN spectrum. We thus apply 
the XSPEC model mekal with absorption {pha(mekal)) for 
an optically-thin thermal plasma with emission lines. Using 
the default abundances (xspec angr), the model provides 
a good fit (x^/d.o.f=29.4/30) with a temperature of 7.3^2 3 
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Figure 6. Continued X-ray count maps and spectra. Upper and lower limits for each count map are 410 and 8.3 counts/arcsec^ for 
G359.933-0.037 (a), 410 and 4.1 counts/arcsec^ for G359. 89-0.08 (b), and 12 and 0.5 counts/arcsec^ for G359. 55+0.16 (c) and G359.43- 
0.14 (d), respectively. 



keV (~8.5l|^ X 10^ K) and NH^IS-Sttl x 10^^ cm-^ The 
best fit normalization of 4.3lg g x 10~^ gives a volume emis- 
sion measure of ~l.lj!o2 cm~® pc^. In § 4.2, we examine 
G359. 942-0. 03 as a ram-pressure confined stellar wind bub- 
ble generated by a massive star. 



3.5 G359.89-0.08 

iLu. Wang, fc Lang I l|2003l ) have previously studied G359.89- 
0.08 (Fig. 6b) in detail from archived Chandra data un- 
der the SNR and PWN scenarios. The SNR case was de- 



termined unlikely in part due to a 10" -30" offset between 
the non-thermal radio emission from the proposed SNR, 
G359.92-0.09, and the X-ray emission from G359. 89-0.08. 
The high X-ray absorption column density also questions 
the SNR scenario as it suggests that G359.89-0.08 lies be- 
hind an enhancement of molecular material, contrary to pre- 
vious work on G359.92-0.09 which indicates that it must 
lie, at least partially, in front of the molecu lar gas (see also 
ILu. Wang, fc Lang il2003l : ICoil fc Holl2000l ). In the previous 
study, they used two observations pointed at Sgr A* (ID 242 
and 1561) with an integrated exposure time of 101 ks. By 
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Figure 8. Region selection for piecewise analyses. The count maps for GO. 13-0. 11 and G359.89-0.08 are overlaid witli regions identifying 
the segments to be used in their respective piecewise analyses. The region selection for GO. 13-0.11 is shown in 8a. The image 8b gives 
the minor axis selection for G359.89-0.08 while 8c gives the major axis selection. 



performing a piecewise analysis of G359.89-0.08, using the 
greater exposure time and counting statistics of this study, 
we may be able to further coustrain the physical model of 
the X-ray emission. We extract spectra from regions along 
the major and minor axes (see Fig. 8b & 8c) and perform 
power law joint fits in which A^h is set as a common pa- 
rameter for the segments. Examining the two fits, wc sec no 
noticeable trend across the minor axis of the feature. The 
outer regions have spectral indexes of 1.3l[5 4 ^'^^ l-^ii i> 
for the SW and NE segments, respectively, while the center 
segment has a spectral index of 1.21q'3, marginally harder 
then either of its neighbors. This is consistent with the PWN 
scenario for synchrotron cooling of the pulsar wind materi- 
als; the outer regions probably represent aging particles in 
a cocoon around a fresh stream from the pulsar. Across the 
major axis, wc sec evidence for possible spectral hardening, 
or flattening, in G359.89-0.08 as shown in Fig. 9 and in a 



contour plot of the feature (Fig. 10), though the overlap 
of the 90% confidence level error bars prevents a definitive 
conclusion. These spectral trends do not provide support- 
ing evidence for the SNR scenario where one would expect 
any detectable systematic trend in the spectral index to be 
along the minor axis while fairly uniform along the major 
axis. Conversely, the apparent observed harderung across the 
major axis is inconsistent with the softening expected for the 
PWN scenario. However, due to the relatively large uncer- 
tainty present in the joint fits and the lack of a confirmed 
pulsar signal, we are unable to judicially rule out the SNR 
case and constrain such spectral evolution. 

3.6 Joint Analysis of Faint Features 

The features GO.223-0.012, G359.55-I-0.16, G359.43-0.14 and 
G359.40-0.08 (Fig. 3a, 6c, 6d and 7, respectively) were ob- 
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Figure 9. Plot of the photon index F across the major axis of 
G359.89-0.08. Positive distance is measured as going northwest 
away from the "head". Error bars are given at 90% confidence. 



-29 05' 45" 




RIGHT ASCENSION (J2000) 



Figure 10. Contour plot of G359.89-0.08. The image, bin size of 
1 pixel (~0.5"), was smoothed using a Gaussian with a FWHM 
of 3 pixels (~1.5"). The upper and lower limits were placed at 
124 and 4.1 counts/arcsec^, respectively. The contour levels were 
placed with a step size of 21 counts/arcscc^ and a starting value 
of 17 counts/arcsec-^. 



served over much lower effective exposure times (in the range 
of RilOO to 150 ks) than those features near the central region 
(ps700 ks to 1 Ms). Due to the lower signal-to-noise of these 
features, the inferred model parameters tend to have rather 
large 90% confidence level error bars (e.g. A^h upper-lower 
bound difference of > 10^"^ and F upper-lower bound differ- 
ence of > 2). To better constrain the parameter bounds, we 
perform joint fits of these 4 'faint' features (faint in terms of 
low signal-to-noise ratios as their fiuxes are still fairly high 
to be seen in the low exposure time). In performing the joint 
fits, we fixed F or A^h to be the same for all features while 
allowing it and the remaining parameters to vary for the 
fit; e.g. fixing A^h to be a common parameter between all 
features and running the fit over the common A'h and each 
feature's F. From these joint fits, we find that the spectral 
parameters of the 'faint' features are consistent over both 
fits. That is, the A^h of the individual features in the F joint 



fit is consistent with the average A^h from the A^h joint fit 
and similarly for F. The spectral parameters for the A'h joint 
fit are given at the bottom of Table 2. Due to the low signal- 
to-noise ratios of the individual features, constraining their 
spectra and morphology is difficult, even with the joint fits. 
Nevertheless, these fits show the relative fiatness present be- 
tween the features and the overall photon index consistent 
with a primarily non-thermal origin. 



4 NATURE OF THE FILAMENTARY X-RAY 
FEATURES 

4.1 Comparison of GC Linear X-ray Threads and 
Known Pulsars 

Outside of the GC, similar linear X-ray features are found di- 
rectly associated with kn own pulsars, providing t he ba sis for 
our PWN interpretation. iKargaltsev fc Pavlov I (|2008l ) have 
collected a number of PWNe detected with Chandra. Many 
of the PWNe listed are affected by their relative motion 
with the ambient medium, e.g. ram-pressure confinement. 
These features closely resemble those in the GC in their mor- 
phological app earance, mos t notably the Mouse PWN (ID 
22 in Kargalts ev fc Pavlov llioos). In deed, a number of the 
features presented in fig. 3 of IKargaltsev fc Pavlov I (|2008| ) 
would appear similar to those found in the GC including 
pulsars J1809-1917, J1509-5850, B1853-H01, G189.23-h2.90, 
G326.12-1.81, a nd G327.15-1.04 (ID 26, 29, 30, 47, 51, 52; 
IKargaltsev fc P avlov 2008) . It should be noted that not all 
X-ray features associated with pulsars can be easily classified 
using a simple model, such as the peculiar case of the pulsar 
B2224+65, also responsible for the Guitar Nebula, and its 
linear X-ray feature which is ~118° offset f rom the pulsar's 
prop er motion (Hui & Becker 2007; Johnso n et al. in prep. I 
I2OO9II . 

Based on known pulsar/X-ray feature associations, we 
are able to derive simple models for many of the GC X-ray 
filaments. While the X-ray counting statistics for many of 
the features studied are indeed low, the power law model 
provides a reasonable fit to the spectra. Out of the 17 fea- 
tures presented, we identify 7 as PWNe based on their 
elongated morphology and spectral parameters. The fea- 
tures GO.13-0.11, G0.03-0.06, G359.97-0.038, G359.97-0.009, 
G359.964-0.052, G359.95-0.04, and G359. 89-0.08 have spec- 
tral indexes in the range of 1.1 to 1.9 with 2-10 keV lumi- 
nosities of 5x10^^ to ^10^"^ ergs s~^, consisten t with typi- 

11200 2: Gacnsl er fc Slane I 
3). The spectral pa- 
rameters of these features support the PWN scenario, 
though the SNR case can not be ruled out, while the low 



cal P WNe (e.g.lGottheff fc Olbert 



I2OO6I : iKaspi. Roberts fc Harding 



signal-to-noise ratios of many of the other features do not 
give enough information to accurately constrain their physi- 
cal models. Moreover, scattering towards the GC makes the 
detection of the radio pulsar signal difficult, preventing a 
definite confirmation of the PWN scenario at present. Detec- 
tion of pulsars in the X-ray band depends on the pulsar's ori- 
entation and mechanism for producing X-ray photons. Ad- 
ditionally, we currently lack the necessary instrumentation 
with enough sensitivity to perform X-ray timing analysis of 
the putative pulsars with the observed faint X-ray fluxes. 
For those features that show emission lines (namely 
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GO.017-0.044 and G359. 942-0.03), the inclusion of a Gaus- 
sian profile with the power-law model reduces the com- 
pared to the power-law only fit with F-test significance 
>98%. It should be noted that, while common practice, 
using the F-test to test for the presen ce of line emis- 
sion hcis a number of short comings (see iProtassov et al.1 
|2002| ). However, we are unaware of a simple, yet statisti- 
cally rigorous, approach that would better determine the 
significance of any lines present. Statistically, thermal mod- 
els, such as Bremsstrahlung, produce satisfactory fits but 
would ha ve unreasonably high tem peratures, i.e. >10 keV 
(see also ILu. Yuan, fc Loul I2OO8I ). The presence of lines 
at either 6.4 keV or 6.7 keV are obvious indicators that 
the features are not PWNe. The 6.4 keV line, in par- 
ticular, is likely due to fluorescent emission from either 
reflection of tr a nsient bright X-ray sources (e.g. Sgr A*; 
iKovama et al. I (Il989l)) or local enhanced low energ y cos- 
mic ray electrons (fyusef-Zadeh. Law fc Wardie1l2002l i. Sim- 
ilar line emissions have been seen in the Galactic Ridge; 
however, the 6.4 keV line emission in the ridge is relatively 
weak, compared to other lines such as the 6.7 keV transi- 
tion, and is consistent with an origin in num erous point-like 
sources (most likely cataclys mic variables, lEbisawa et al. I 
l2008l : iRevnivtsev et al.1l2006l ). 



4.2 G359. 942-0. 03 as a Ram-Pressure Confined 
Stellar Wind Bubble 

We find a point-like near-IR source that may be asso- 
ciated with G359. 942-0.03 (Fig. 5c) in the 2MASS cata- 
log (ID 17453582-2900050, ISkrutskie et al. I [2OO6 ) and in 
a recent HST/NICMOS s urvey of the GC at 1.90 ^im 
(|Wang et al. in prep. II2009I ). The 2MASS counterpart has 
an apparent H-band magnitude of 10.85±0.06 and an H- 
K color of 1.67 mag, in agreement with those observed for 
massive stars in the GC region. The 2MASS source has 
an angular separation of .71" with respect to the approx- 
imate center of the "head" of G359. 942-0.03. Within 30" of 
G359.942-0.03, there are a total of 48 2MASS sources re- 
lating to a P-statistic (chance of random association, calcu- 
lated using P=l-e where n is the 2MASS source den- 

sity and 6 is the dist a nce b e tween potent i al cou nterparts, 
see also iBorvs et al.1 (|2004l ): iPope et al. I (|2006l )) of ~3% 
for the potential 2MASS counterpart. If one considers the 
size of the extended X-ray emission, then this P-statistic in- 
creases to ~50%; however, if this model is correct, we should 
only consider the region where the massive star would be lo- 
cated; e.g. in the "head", leading to a robust identification 
for the 2MASS counterpart. Interestingly, the color excess 
E(H-K) = 1.72 (assuming (il-K) int r insic ~-0.05 for early-type 
stars (see also lFiger et al. I ll999al : [Panagia1ll973t )') is signif- 
icantly lower than that inferred from the X-ray absorption. 
This discrepancy is expected if the star has a strong wind, 
such Eis a Wolf-Rayet star, which can produce a free-free 
emission enhancem ent along with H/He line s in th e H and 
K bands (see also iFiger et al. 1 1 19971 . Il999al . I2OO2I '). As an 
example , the brightest mas sive star identified in the Arches 
Cluster l|Figer et al. II2002I ) (2MASS ID 17455043-2849215) 
has apparent H and K band magnitudes of 9.5 and 8.6, 
respectively, giving E(H-K)=0.95. Additional redding may 
come from the large variations in extinction around the 



GC over large distances. If we can then assume that this 
source is a massive star with strong winds, one may ex- 
pect to see an extend ed emission region in the Paa map of 
IWang et al. in prep. ] (12009 ): however, the expected extinc- 
tion based on the X-ray column density would place the Paa 
flux well below the detection limit of the survey. This indi- 
cates that while the presence of strong winds associated with 
a massive star can not be confirmed, this scenario remains 
plausible for the 2MASS counterpart of G359. 942-0.03. 

This ram-pressure confined wind bubble model is driven 
by the presence of the massive star as well as the cometary 
morphology and the apparent thermal spectrum of the X- 
ray emission (§ 3.4). In Fig. 5c, we can see fairly well con- 
strained extended X-ray emission, with apparent dimensions 
of ~2"x6" or .08X.23 pc, with a bright "head" region, con- 
sistent with our definition of a cometary feature. In this 
model, the motion of the star creates a bow shock ahead of it 
that sweeps up material from the ISM, sim ilar to the model 
presented bv lvan Buren fc McCravl (|l988t ). We may infer a 
number of model parameters for this scenario based on the 
properties of the X-ray-emitting gas (§ 3.4). First, the strong 
stellar wind from the massive star is expected to be heated in 
a reverse shock to a temperature of T ~ 6 x lO'^ K(ww/2)'^, 
where v„ is the stellar wind velocity (in units of 10^ km 
s~^). The possible range of stellar wind velocities under 
the massive star scenario, typically of 1-2x10^ km s"'^, can 
then produce temperatures high enough to match the value 
(^8.5^2 g ^ lO'^ K) inferred from the spectral analysis. Sec- 
ond, the measured volume emission measure (~l.ll'2 cm~^ 
pc^) provides an estimate on the density of the shocked wind 
material of rie ^ 30 cm"'', assuming a cylindrical volume of 
the observed width and length (Table 1). Third, constraints 
may be placed on the stellar velocity from balancing the 
ram and shocked wind pressures, 2nekT ~ 10^'' fiHrrinnov^ , 
where is the stellar velocity (km s~^), /lu is the average 
molecular weight per hydrogen atom (^h=0.6 for approxi- 
mately solar metallicity) and no is the density of the ambient 
ISM (cm""'). This balance of the ram-pressure also serves to 
determine the size of the bow shock region. Van Buren fc 
McCray (1988) showed that the distance to the contact dis- 
continuity of the bow shock is ~ 1.5Zi where h is the length 
scale characterizing the point where the ISM ram pressure 
equals that of the stellar wind. The length scale h is then 
given by 



h = (56 pc)M'/"u 



1/2 1/2 

/ .1, ' T 



1/2 -1/2 



(1) 



whe re M is the mass loss rate in units of 10 ^ Mq yr ^ (see 
also lvan Buren fc McCrav"1ll988l . eq. 1). Since the measured 
width of G359.942-0.03, ~0.08 pc taken across the "head" 
region, should be ^ 3/i, as the shocked wind gas is primarily 
responsible for the X-ray emission, we have a stellar velocity 
estimate of 



(27 km s-i)Af^/2«y2(no/10^)-^/^ 



(2) 



From the above equations, we then infer no ~ 1.5 — 
7 X 10^ assuming M = 1 and Hw = 2 or ~ 39 - 84 km s~^ 
for no « 10* cm~^, consistent with what may be expected 
from a run-away massive star into a dense cloud. Based on 
typical parameters for massive stars, this ram-pressure con- 
fined wind bubble model is capable of explaining the thermal 
spectrum, luminosity and well collimated, comet-like mor- 
phology of the X-ray feature G359. 942-0. 03 along with its 
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apparent association with a massive star. We should caution 
readers that this conclusion, while reasonable, is not definite 
due to the number of assumed parameters. Further analy- 
sis, including deeper infrared imaging and spectroscopic ob- 
servations, can aid in determining the assumed values and 
provide confirmation (or contradiction) to this ram-pressure 
confined stellar wind bubble model. 



5 SUMMARY AND FINAL REMARKS 

To conclude. X-ray features in the GC likely have a hetero- 
geneous origin. Non-thermal features are typically produced 
through synchrotron emission of PWNe or magnetohydro- 
dynamical shock fronts from SNRs. For the PWN model, one 
would expect elongation due to either the pulsar movement 
in the ISM or magnetic field confined fiows of the pulsar 
wind materials. In such a case, one would expect to see spec- 
tral softening across the direction of elongation along with 
typical photon indexes of r=l.l-2.4. For the SNR case, the 
direction of elongation would be caused by shock fronts and 
any measurable spectral change would be expected across 
the width of the feature. Of the features presented, 7 exhibit 
characteristics of PWN, although the SNR case has not been 
completely ruled out. Their non-thermal spectra as well as 
typical comet-like morphology are consistent with the ram- 
pressure confined PWN model. The lack of any detected 
pulsar signal from these PWNe prevents a definitive conclu- 
sion on their origin. Some of these pulsars may be detected 
in future pulsar searches using high frequency and high res- 
olution radio observations. In the case of G359. 942-0. 03, the 
model of a ram-pressure confined stellar wind bubble best 
explains the detection of the 6.7 keV He-like Fe emission and 
steep spectral shape. The near-IR images show an apparent 
counterpart of G359. 942-0. 03 consistent with a massive star 
where the lack of extended near-IR emission could be due 
to a combination of extinction and inefficiency in ionizing 
photon conversion. 

In addition to their physical models, the X-ray features 
provide information into the underlying structure of the GC. 
Featur es such as GO. 017 - 0. 044 and the "clumps" as stud- 
ied bv ILu. Yuan, fc Lou I (|2008l ) map out the history of the 
GC (e.g. the possible radiative illumination of GO. 017-0. 044 
from Sgr A*). The orientations of PWN-like features help 
to map the magnetic, ga s and stellar dynamics of the GC. 
ILu. Yuan, fc Loul l|2008h predicted that the orientation of 
comet like features is due to a strong galactic nuclear wind 
comparable to typical pulsar velocities oriented away from 
Sgr A*. Of the features present within a 20'x20' image cen- 
tered on Sgr A* (see Fig. 2), this effect can be seen in 5 
out of the 7 PWN candidates (excluding G359. 95-0.04 and 
G359.89-0.08) and G359.933-0.037. Those features outside 
the 20'x20' region along with the features G359. 964-0.03, 
G359.95-0.04 and G359.89-0.08 do not necessarily follow this 
claim; though this could be explained by variations in the 
vector field of the galactic wind as a result of regions of vary- 
ing density. Additionally, the filament-like features could be 
regulated by local ordered magnetic field lines as with G0.13- 
0.11's morphology. These features can thus be used in more 
in-depth studies of GC dynamics, beyond the quick analy- 
sis just presented, in combination with deep observations at 



other wavelengths; e.g. radio polarization and magnetic field 
tracers. 

For the GC observations, few X-ray filaments were iden- 
tified in regions outside the central 20'x20' region. These re- 
gions have integrated exposure times significantly less than 
the central region. In addition, few of the identified features 
within the central 20' x 20' region were found in initial obser- 
vations with exposure times <200 ks. Thus, those presented 
in this study may represent only a tip of the "iceberg" of 
X-ray features. Indeed, some may be aligned such that they 
appear like point sources with their direction of elongation 
along the line of sight. Based on the ave rage star formatio n 
rate of the GC (^ W2!_Mq yr'^ pc~ MFiger et al. II2004I ). 
ILu. Yuan, fc Loul (|2008t ) predicted ~40 pulsars within 7 pc 
of the GC with ages < 3 x 10"" yr assuming an average stel- 
lar mass of IOMq. This is roughly consistent with our 13 
PWN candidates if we consider that projection effects will 
limit us to ~71% of PWN, assuming the extended X-ray 
emission of a PWN is identified if the angle between the ex- 
tended emission and the line of sight lies between 7r/4 and 
37r/4, along with factors including sensitivity limits and pul- 
sar population statistics. With higher sensitivity and spatial 
resolution in multiple wavelength bands, including X-ray, ra- 
dio and near-IR, the identities of these X-ray features can 
be more accurately determined. 
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Table 1. Morphological Properties of Identified X-ray Features 



ID 


Shape 




Size 

pc 




Orientation 


Notes 


GO. 223-0. 012 


fil cimcnt jiry 


(] 


23 X 3 


12 


ESE-WNW 


Faint 


GO. 13-0. 11 


comctary 





28x1 


29 


ESE-WNW 


Curved feature 


GO. 03-0. 06 


comctary 





27x1 


IB 


SE-NW 


Curved 'tail' 


GO. 017-0. 044 


filamentary 





08x0 


50 


ESE-WNW 


Faint 


GO. 007-0.014 


filamentary 





12x0 


39 


ENE-WSW 


Faint 


G359. 97-0.038 


cometary 





23X0 


54 


SW-NE 




G359. 97-0. 009 


cometary 





08X0 


39 


SSE-NNW 




G359. 964-0. 052 


filamentary 





08X0 


50 


NNE-SSW 




G359. 96-0. 028 


filamentary 





12x0 


39 


SE-NW 




G359.95-0.04 


comctary 





08x0 


31 


NNE-SSW 




G359. 942-0. 03 


cometary 





08x0 


23 


ENE-WSW 


Faint 


G359. 94-0.05 


filamentary 





08x0 


39 


ESE-NWN 




G359. 933-0.037 


cometary 





08x0 


27 


ENE-NWN 


Faint 


G369. 89-0.08 


cometary 





31X0 


93 


SE-NW 




G359. 55-1-0.16 


filamentary 





31X2 


18 


E-W 


Faint 


G359.43-0.14 


romctarv 


1) 


l.'')X() 


83 


N-S 


Faint 


G:i.")!J.4u-().0.'i 




(J 


20X1 


U7 


S-N 


Faint.. (. 'ui \ fd Iraturc 



Note. — Sizes are taken from apparent surface brightness as compared to 
the local background. Orientation is based on North and East defined as 
up and left in Chandra images, respectively. Notation for orientation then 
follows that of the cardinal directions. "Faint" refers to features that have 
low signal-to-noisc or low surface brightness when compared to the respective 
local background. 



Table 2. X-ray Spectral Properties of Identified Features 



ID 




r 






EW 


Line Energy 






(10^2 cm-2) 




(10~^* ergs/s/cm^) 


(10^^ ergs/s 


(keV) 


(keV) 




GO. 223-0.012 


_ ,4-79.8 
°-^-5.1 




15 


15 






6.7/8 


GO. 13-0.11 




l.52l°;| 


34 


40 






96.0/100 


GO. 03-0.06 


6.3l-_«, 


1 1+°^* 

-^■-^-0.4 


12 


13 






118.9/119 


GO. 017-0.044 




-0 7+'^-° 
"■'-0.6 


4 


3 




6 4+°-^ 
"■*-o.i 


15.4/16 


GO. 007-0.014 


^■'-5.7 


1 0+3-* 


2 


2 






24.3/20 


G369. 97-0.038 


11 7+^-^ 
' -1.9 


1 4+0-8 
-^-■^-0.3 


11 


15 






109.7/137 


G369. 97-0.009 


n fj-t-8.8 


1 2+°-'' 
'-■''-0.6 


4 


5 






52.5/59 


G359. 964-0.052 


11 1 + ^^* 


'^■"-0.3 


23 


37 
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Note. — Spectral parameters obtained through XSPEC model fitting (pha(po) or pha(po+Ga)) to each of the iden- 
tified X-ray features. In order, the parameters are the hydrogen gas column density A^h, X-ray photon index F, 
observed X-ray flux and unabsorbed X-ray luminosity in the 2.0-10.0 keV band Fx and Lx, equivalent width and 
energy centroid of the einission line if present, and pi^r degree of freedom. Luminosities are based on the assumed 
distance of 8 kpc to the Galactic center. Spectral parameters for "faint" features are based on the joint fit model 
with ATh as a common parameter. Errors are given to the 90% confidence interval. 



